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Abstract 

We explore the collective electromagnetic response in atomic clusters of various sizes and ge- 
ometries. Our aim is to understand, and hence to control, their dielectric response, based on a 
fully quantum-mechanical description which captures accurately their relevant collective modes. 
The electronic energy levels and wave functions, calculated within the tight-binding model, are 
used to determine the non-local dielectric response function. It is found that the system shape, 
the electron filling and the driving frequency of the external electric field strongly control the res- 
onance properties of the collective excitations in the frequency and spatial domains. Furthermore, 
it is shown that one can design spatially localized collective excitations by properly tailoring the 
nanostructure geometry. 
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I. INTRODUCTION 



Recent advances in nanoscience have created a vast number of experimentally accessible 
ways to configure atomic and molecular clusters into different geometries with strongly vary- 
ing physical properties. Specifically, exquisite control of the shape and size of atomic and 
molecular clusters has made it now possible to investigate the collective electromagnetic re- 
sponse of ultra-small metal and semiconductor particles.- 1 ^ The aim of this study is to model 
and examine plasmonic excitations in such structures, and thus to gain an understanding of 
the quantum-to-classical crossover of collective modes with increasing cluster size. There is 
obvious technological relevance to tunable collective modes in nanostructures. For example, 
surface plasmon resonances in metallic nanospheres and films have been found to be highly 
sensitive to nearby microscopic objects, and hence are currently investigated for potential 
sensing applications.- In this context, it is desirable to design customized nanostructures 
with specifically tailored resonance properties,- and this study is intended to be a first step 
into this direction. 

It is natural to expect that in many cases the electromagnetic response of nanoclusters 
is considerably different from the bulk. In particular for very small clusters, the quantum 
properties of electrons confined in the structure need to be taken into account. Moreover, 
unlike in the bulk, the coupling between single-particle excitations and collective modes can 
be very strongly affected by its system parameters. This exponential sensitivity opens up 
excellent opportunities to optimize the dielectric response via tuning the cluster geometry 
and its electron filling. For example, by proper arrangement of atoms on a surface one can 
design nanostructures with controllable resonances in the near-infrared or visible frequency 
range.- A possible application of such nanostructures is the creation of metamaterials with 
negative refractive index at a given frequency. Furthermore, since geometry optimization 
of bulk resonators has demonstrated minimization of losses in metamaterials 5 , it is also 
interesting to investigate the effect of the nanostructure shape on the loss function at a 
given resonance frequency. 

To approach this problem, in this study we investigate the formation of resonances in 
generic systems of finite conducting clusters, and examine how their frequency and spatial 
dielectric response depends on the system size and geometry. In particular, the non-locality 
of the dielectric response function in these structures is important, and will therefore be 
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properly accounted for. A similar analysis for the case of small metallic nanostructures was 
performed recently, using an effective mass approximation.- Here we focus on the opposite 
limit, namely we assume that electrons in the cluster can be effectively described using a 
tight-binding model.- Because of the more localized nature of the electronic wave functions 
in this model the overall magnitude of the collective modes are expected to be strongly 
suppressed as compared to metallic clusters. 

This paper is organized in the following way. In the next section we introduce the model 
and method. For a more extended discussion, the reader is referred to RefA Results for 
the induced energy as a function of the driving frequency of an externally applied electric 
field and the corresponding spatial modulations of the charge density distribution func- 
tion are discussed in the following section. Finally, a discussion of possible extentions and 
applications is given in the conclusion section. 

II. MODEL 

The interaction of electromagnetic radiation with nanoscale conducting clusters is con- 
ventionally described by semi-classical Mie theory- This is a local, continuum-field model 
which uses empirical values of the linear optical response of the corresponding bulk mate- 
rial, and has been applied in nanoparticles to describe plasmon resonances.- However, such 
a semi-empirical continuum description breaks down beyond a certain degree of roughness 
introduced by atomic length scales, and thus cannot be used to describe ultra-small systems. 
In addition, near-field applications, such as surface-enhanced Raman scattering,— are most 
naturally described using a real-space theory which includes the non-local electronic response 
of inhomogeneous structures. Therefore, we will use a recently developed self-consistent and 
fully quantum-mechanical model which fully accounts for the non-locality of the dielectric 
response function.- 

Specifically, to identify the plasmonic modes in small clusters we calculate the total 
induced energy due to an applied external electric field with driving frequency u, and scan 
for the resonance peaks. The induced energy is determined within the non-local linear 
response approximation. 

To keep the computational complexity of this procedure at a minimum, we use a one- 
band tight-binding model to obtain the electronic energy levels Ei and wave functions ^(r) 
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linear combination of of s-orbitals 



i>i( T ) = ^ViMr-Rj), (1) 



h3 



where (p(r — Rj) is the wave function of a s-orbital around an atom localized at position Rj 
and otij are the coeficients of the eigenvector (with energy Ej) of the Hamiltonian, which 
has the matrix elements 



fi i = j 



(<p(r-B i )\H\<p{T-B i )) = { 



—t i,j n.n. (2) 
otherwise 



Here t is the tight-binding hopping parameter. The Hamiltonian matrix is diagonalized 
using the Householder method to first obtain a tridiagonal matrix and then a QL algorithm 
for the final eigenvectors and eigenvalues. 11 

Once the electronic wave functions have been obtained, it is possible to calculate the 
dielectric susceptibility %(r,r', oS) via 

X (r,r',u) = ^^^- <(r)^(rO^(rO^(r). (3) 

The induced charge density distribution function is then obtained by 

Pind(r,uj) = J x(r,r',u;)(0i nd (r / ,a;) + <p ext (r' ,u))dr' , (4) 

where in turn the induced potential is given by 

4> ind [r,u)= / — —dr . (5) 

J \ r — r | 

We avoid the large memory requirement to store x( r ) r 'i ^) by calculating the induced charge 
density distribution iteratively via 

p md (v^) = J {E i~ fiEi) . ^(r)^-(r) / ^(r')^(r',u;)^(r')dr', (6) 
^ E, t - Ej - lo - ij J J 



with (j) tot {r',u) = c/)i n d(r' , uj) + (j) ext {r' , u) . The integrals are evaluated using a 4* h order 
formula obtained from a combination of Simpson's Rule and Simpson's 3/8 Rule. Eqs. [5] and 
[6] are solved self-consistently by iterating (j) ind (r,u) and p ind (r,uj). This procedure typically 
converges in 3-8 steps when starting with (f) ind (r,u) = 0, depending on the proximity to 



a resonance and on the value of the damping constant 7, which throughout this paper is 
chosen as 7 = 0.081 A much better performance can be achieved when the initial (pi n d(r,uj) 
is taken as the solution of a previously solved nearby frequency. Upon its convergence, the 
frequency and spatial dependence of the induced electric field and the induced energy are 
obtained using 



The observed resonances in the induced energy and charge density distribution at certain 
driving frequencies of the applied electric field correspond to collective modes of the cluster. 
In the following, the local induced charge density distribution is used for analyzing the 
characteristic spatial modulation of a given plasmonic resonance. 

III. RESULTS 

Let us first focus on the dielectric response function in linear chains of atoms, with the 
intent to identify the basic features of their collective excitations. The frequency dependence 
of the induced energy in such systems, exposed to a driving electric field along the chain 
direction, is shown in FigJT](a). It exhibits a series of resonances, which increase in number 
for chains with increasing length. As observed in the spatial charge density distribution, e.g. 
shown for the 5-atom chain in Figfjjb), the lowest peak corresponds to a dipole resonance. 
When increasing the system size, the dipole peak moves to lower frequencies, which is 
the expected finite-size scaling behavior. The resonances at higher frequency correspond 
to higher harmonic charge density distributions. For example, in FigJT](c), we show the 
charge density distribution corresponding to the highest frequency resonance of the 6-atom 
chain. In contrast to the dipole resonance, these modes show a rapidly oscillating charge 
density distribution, and thus have the potential to provide spatial localization of collective 
excitations in more sophisticated structures. While an extension to much larger chains 
is numerically prohibitive within the current method, the finite-size scaling of the observed 
dielectric response of these clusters is consistent with the ID bulk expectation of a dominant 
low-energy plasmon mode, coexisting with a particle-hole continuum of much weaker spectral 
intensity. 
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FIG. 1: Longitudinal modes in atomic chains, (a) Decimal logarithm of the total induced energy 
(artificially offset) as a function of the frequency of an external electric field which is applied along 
the direction of the chain. The resonance peaks correspond to different modes, (b) Induced charge 
density distribution for the lowest energy mode at uj = 0.73t in the 5-atom chain, (c) Induced 
charge density distribution for the highest-energy mode at uj = 3.56t in the 6-atom chain. 

In order to study the transverse collective modes we apply an external electric field 
perpendicular to ladder structures made of coupled linear chains of atoms.-^ Figj2](a) shows 
that for every chain size there are two resonance peaks for the total induced energy, the higher 
energy is an end mode, as shown in Figsf2](b) and (c) for the 3 and 5-atom double chains 
respectively, whereas the lower energy peak corresponds to a central mode, as displayed in 
FigfS^d) for the 6-atom double chain. It is also confirmed that as the length of the chain is 
increased, the central mode gets stronger relatively to the end mode, which is the expected 
behavior for bulk vs. surface excitations. These results are in agreement with the findings 
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FIG. 2: Transverse Modes in coupled chain structures, (a) Logarithm of the total induced energy 
(artificially offset) as a function of the frequency of an external electric field applied transversely to 
the chain. The low energy mode is central, the analog of a bulk plasmon, and the high energy mode 
is located at the surface, the analog of a surface plasmon. (b) Induced charge density distribution 
for the mode at u = 4.91t in the 3-atom double chain, (c) Induced charge density distribution for 
uj = 5.41i in the 5-atom double chain, (d) Induced charge density distribution for u = 1.93i in the 
6-atom double chain. 

of Ref A 

Let us next examine what happens when the direction of the external electric field is 
varied. Figj3] shows the dielectric response of a 4 x 6-atom rectangular structure for different 
angles incidence directions of the applied field. When the field is parallel to one of the edges 
(6 = 0°or 9 = 90°), the response is essentially that of a single chain with the same length, 
shown in Fig{T](a). Also the induced spatial charge density modulations are analogous to 
those of the correspondent linear chain, which can be seen in Figs. Gl^b) and (c). At 
intermediate angles the response is a superposition of the two above cases, changing gradually 
from one extremum to the other as the angle is changed. Notice for instance that as the 
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FIG. 3: Dependence on the direction of the applied electric field, (a) Logarithm of the total 
induced energy (artificially offset) as a function of frequency of external electric fields applied to 
a 4x6-rectangle at different incident angles. 9 = 0°when the field is parallel to the 4-atoms edge, 
and 9 = 90°when it is parallel to the 6-atoms edge, (b) Induced charge distribution for 6 = 0°and 
<jj = 3.15£. (c) Induced charge density distribution for 9 = 90°and u = 2.15t. 

angle increases, the peak at the same frequency of the 4-atoms dipole resonance diminishes, 
while simultaneously another resonance is formed at the frequency of the dipole mode of a 
6-atoms chain when the angle is tuned from 9 = 0° to 9 = 90°. For 9 = 0° there is only the 
peak at the frequency of the 4-atom chain dipole resonance, whereas for 9 = 90° only the 
dipole peak corresponding to the 6-atoms dipole frequency is present. The superposition 
of the responses from each direction is a consequence of the linear response approximation 
employed, since the response is a linear combination of those obtained from each direction 
component of the external field. 

Next, let us analyze the dependence of the resonance modes on the number of electrons 
in the cluster. FigJUa) shows significant changes in the response of a 9-atom chain with the 
external field applied along its direction. In particular, it is observed that the response is 
stronger when there are more electrons in the sample, a quite obvious fact since there are 
more particles contributing to the collective response. Moreover the resonance frequencies 
of lower modes increase with the number of electrons, which can be understood as a con- 
sequence of the one-dimensional tight-binding density of states being smallest at the center 
of the band. Hence the energy levels around the Fermi energy are more sparse in the finite 
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FIG. 4: Variation of the number of electrons, (a) Logarithm of the total induced energy as a 
function of the external electric field frequency. The number of electrons N e i in a 9-atom chain 
is varied, (b) Induced charge density distribution for N e i = 1 at u = O.SQt. (c) Induced charge 
density distribution for N e i = 9 at uj = 0.53i. 

system, and therefore the excitations require larger frequencies at half-filling. The same 
does not hold for higher frequency modes since these correspond to transitions between the 
lowest and highest levels for any number of electrons in the sample. Therefore these modes 
have the same frequency, independent of the electronic filling. Higher filling also allows 
the induced charge density to concentrate closer to the boundaries of the structure, as a 
comparison between FigsJH(b) and (c) demonstrates. FigJU^b) shows that a 9-atoms chain 
with N e i = 1 electron has its induced charge density localized around the center of the 
chain. In contrast, FigJH(c) displays the induced charge density localized at the boundaries 
of the same structure with N e i = 9. This concentration closer to the surface happens be- 
cause higher energy states have a stronger charge density modulation than the lower energy 
ones. Therefore the induced charge density is more localized for higher fillings, because at 
low fillings the excitations responsible for the induced charge density are between the more 
homogeneous lower energy levels. This can be interpreted as a finite-size rendition of the 
fact that by increasing the electronic filling one obtains the classical response with all the 
induced charge density on the surface of the object. 

Access to high energy states is very important for achieving spatial localization of the 
induced charge density, as the next example shows. In order to find a structure with spatially 
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FIG. 5: Connect between two chains, N e i = 1. (a) Logarithm of the total induced energy as a 
function of the frequency of an external electric field applied to two 8-atom chains with an extra 
atom connecting them at the center, (b) Induced charge density distribution for uj = 2.28t. (c) 
Induced charge density distribution for u> = 4.36t. 

localized plasmons we consider two parallel 8-atoms chains connected to each other by an 
extra atom at the center. When an external electric field is applied transversely to the 
chains, the electrons are stimulated to hop between them, but this is only realizable through 
the connect, therefore the plasmonic excitation is sharply localized around it. Figj5](a) shows 
the response of two 8-atoms chains, Figj5]^b) and (c) show respectively the induced charge 
density for the dipole and the highest modes. It is seen that the lowest frequency mode has 
a much longer spread along the chains, whereas the high frequency plasmon is sharper, since 
it corresponds to excitations to the highest energy state that has a large charge ondulation 
as it was pointed out before. 

Let us finish by analyzing the dependence of the various dielectric response modes on the 
distance between atoms in the chain. The dipole moment of the chain is proportional to 
its length, and consequently also proportional to the distance between neighboring atoms. 
Hence one can expect that the strength of the dielectric response is also proportional to 
the atomic spacing. However, higher frequency modes require a larger mobility of electrons 
since they need to move quickly along the chain in order to produce the charge oscillations 
of the mode. In other words, the charge- density waves can only have frequencies at most 
equal to the wavering of the charge carriers. Hence the high frequency modes are suppressed 
for systems where electrons lack the mobility necessary for the execution of their oscillation. 
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FIG. 6: Variation of the distance between neighbor atoms a. Logarithm of the total induced energy 
as a function of the frequency of an external electric field applied to 7-atom chains with different 
spacing between atoms a in units of the Bohr radius r#. The frequency unit is £3 the tight-binding 
hopping parameter for a = 3rs- 

In the tight-binding model, the electronic mobility stems from the tunneling of electrons 
from one atom to another. The tunneling rate depends on the overlap of the atomic orbitals 
located in different sites, and moreover this overlap decays exponentially with the spacing 
between atoms. Therefore the high frequency modes are weaker for long spaced chains, 
while the slow modes get stronger. This fact is demonstrated in FigE] where the response 
of 7-atom chains with different atomic spacings are shown. The tight-binding hopping 
parameter t changes with the atomic spacing a, and here we considered a generic^ power- 
law dependence t ~ a~ 3 . The figure clearly shows that the strength of the slowest mode 
is reduced for shorter spacings. The opposite is true for the two fastest modes, while the 
intermediate modes have nearly no change. 

IV. CONCLUSION 

In conclusion, we have analyzed the evolution of plasmonic resonances in small clusters as 
a function of the system shape, applied external fields, electron filling and atomic separation. 
Using a fully quantum-mechanical, non-local response theory, we observe that longitudinal 
and transverse modes are very sensitive to these system parameters. This is reflected in 
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their frequency, oscillator strength and the spatial modulation of the induced charge density. 
Specifically, we identify bulk and surface plasmonic excitations which can be controlled in 
amplitude and frequency by the cluster size. Furthermore, we observe a non-trivial filling 
dependence, which critically depends on the electronic level spacing in a given structure. 
We also find that changes in atomic spacings strongly affect the electron mobility in these 
structures, with a very different impact on low-energy vs. high-energy modes. And we 
see that changing the position of a single atom in a nanostructure can completely alter its 
collective dielectric response. This strong sensitivity to small changes is the key to controlling 
the modes of ultra-small structures, and it can thus become the gateway to a new generation 
of quantum devices which effectively utilize quantum physics for new functionalities. 
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